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We describe a general method to study the ground state phase diagram of electronic models on 
chains whose extended Hubbard hamiltonian is formed by a generalized permutator plus a band- 
controlling term. The method, based on the appropriate interpretation of Sutherland's species, yields 
under described conditions a reduction of the effective Hilbert space. In particular, we derive the 
phase diagrams of two new models; the first one exhibits a band-controlled insulator-superconductor 
transition at half-filling for the unusually high value U c = 6t; the second one is characterized by a 
filling-controlled metal-insulator transition between two finite regions of the diagram. 
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Metal-Insulator and Insulator-Superconductor Transi- 
tions in chain systems of interacting electrons have re- 
cently become a matter of great interest for the physics 
of new compounds and devices Q. Although many ex- 
perimental data are nowadays at disposal, an important 
open question of this issue is to determine the Hamilto- 
nian that could fairly describe these kinds of transitions. 
The task is particularly difficult just owing to the low di- 
mensionality, which causes usual mean-field and pertur- 
bative approaches to often fail in providing reliable pre- 
dictions. Fortunately, 1-dimensionality allows to exploit 
exact analysis techniques which can provide -although 
only for some particular cases- rigorous informations on 
the structure of the ground state and on low-energy ex- 
citations. For this reason a probative test for theoretical 
models is the comparison between experimental results 
and theoretical predictions on the ground state phase di- 
agram. 

The ground state is usually given as a function of the fill- 
ing n, i.e. the number of effective carriers, and of a 'band- 
parameter', which indicates the intrinsic unit of energy of 
the system (its actual definition depends on the theoreti- 
cal approach envisaged, see below). One can thus distin- 
guish between filling-controlled (FC) and band-controlled 
(BC) transitions, according to which kind of parameter 
(number of carriers or energy scale respectively) is tuned 
to let the transition occur. Both kinds of transitions are 
very important in practical applications: BC transitions 
are relevant for instance in Vanadium oxides, where one 
can modify the band-width through hydrostatic pressure 
on the sample; FC transitions are frequent in perovskite- 
like materials such as Ri^^A^TiOa (R=rare-earth ion, 
and A=alkaline-earth-ion), as well as in hole-doped com- 
pounds like Lai_ x Sr :E Cu02.5. In order to describe these 
materials, at least as far as their low energy excitations 
are concerned, where single band picture are often reli- 



able, the class of extended Hubbard models provides an 
interesting starting point. For these models, which in- 
volve strong electronic correlations, the band-parameter 
is usually [El taken to be the on-site Coulomb repulsion 
(U), instead of w = At (t = hopping amplitude), the lat- 
ter being typical in mean-field approaches. A number of 
exact results have been obtained in terms of n and U. For 
the ordinary Hubbard model, a FC metal-insulator-metal 
transition has been shown to occur at half filling (n = 1) 
for any U > (U being the on-site Coulomb repulsion); 
on the contrary, no BC transition takes place for U > 0. 
More recently, some models ( || and Q]) were solved in 
which a BC Insulator-Superconductor Transition occurs 
at half- filling at finite values of U > 0, while the usual FC 
Metal-Insulator-Metal transition takes place for n = 1 
and U > U c . 

At the best of the authors' knowledge, no detailed investi- 
gation has been devoted to either of the following issues: 
for the BC transitions it has not been pointed out yet 
what interaction terms are relevant to tune the critical 
value U c at which the transition occurs: this is quite im- 
portant because U c can assume different values according 
to the chemical structure of the material. Secondly, for 
the FC transitions, all the above models provide an in- 
sulating state only at half filling; on the contrary, doped 
materials exhibit an insulating phase for a finite region of 
filling values. In this letter we examine the above subjects 
providing the exact ground state phase diagram of some 
1-D electronic models. In particular, we obtain rigorous 
results which allow us to both discuss the dependence 
of U c of BC transitions on the Hamiltonian parameters, 
and to find a FC metal-insulator transition between two 
finite regions of the phase diagram. 

We consider here the most general 1-band extended 
isotropic Hubbard model preserving the total spin and 
number N of electrons, which reads 
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In (|l]) Cj a ,Cj.cr are fermionic creation and annihilation 
operators on a 1-dimensional chain with L sites, a G {| 
,1}, fij^ = c^ a c jjCr , fij = X) CT ^>. and (j, fc) stands for 
neighboring sites, t represents the hopping energy of the 
electrons (henceforth we set t = 1), while the subsequent 
terms describe their Coulomb interaction energy in a nar- 
row band approximation: U parametrizes the on-site re- 
pulsion, V the neighboring site charge interaction, X the 
bond-charge interaction, W the exchange term, and Y 
the pair-hopping term. Moreover, additional many-body 
coupling terms have been included in agreement with |3| : 
X correlates hopping with on-site occupation number, 
and P and Q describe three- and four-electron interac- 
tions. In the following we shall identify the 4 physical 
states | |), | J,}, |0) and | if) at each lattice site with 
the canonical basis e Q of C 4 , and denote nt = N^/L; 
ni = Ni/L; n Q — N a /L; n-j-j. = N^JL the densities of 
these 4 species of physical states. 

In H it has been shown that, by fixing all the coupling 
constants of (Q) to appropriate values, one can rewrite 
TL as a generalized permutator (GP) between neighbor- 
ing sites (minus some constant terms). Here we add to 
the latter a further arbitrary term UYlj nj^fiji, which is 
easily proved to commute with the GP. In matrix repre- 
sentation, the Hamiltonian ([!]) that we consider reads 

H = — Tl a p + UNfi — const.terms (2) 

Q>/3 

where U a p acts as a GP II (see below) whenever 2 neigh- 
boring sites of the chain are occupied by e a and ep, oth- 
erwise it gives zero. The constant terms are of the form 
UN-ii + p,N + cl, where U, ft and c are fixed values. 
The purpose of this letter is to show how to investigate 
the ground state phase diagram of (|2|) as a function of 
the band parameter U and the filling of the carriers n. 

Let us first recall some basic properties of the GP's. 
With respect to an ordinary permutator, a generalized 
permutator can either permute or leave unchanged the 
states of the 2 neighboring sites (including a possible ad- 
ditional sign); explicitly: 

LT (e a ® e/9 ) = (e Q <g> e ) + 6° a/3 {e <g> e a ) (3) 

where and 9° a ^ are two discrete valued (0, -1 or 1) 
functions determining on II the positions and the signs of 
the diagonal and off-diagonal entries respectively. Also, 
9 d and 9° are 'complementary', i.e. \9^\ = 1 - \9° a ^\, 
so that II has only one non-vanishing entry for row or 
column. Moreover, 9° a p = 0^ a , hence LT is a symmetric 
matrix. The set of couples of subscripts (a, (3) for which 
Kp + ( res P-^ + °) is denoted by A d (resp.-4°). It is 
easily seen that A d is always of the form A d = Uj Sj x Si, 
the 5i's being disjoint subsets of the set S = {1,2,3,4}. 
By varying the functions 9 d and 9° one obtains different 
kinds of GP's. 



In order to determine the ground state phase diagram 
of (0), the difficult task is to calculate the contribution to 
the energy of the first term. To solve this issue, one can 
reconduct it to a problem defined in a smaller Hilbcrt 
space. Indeed it can be seen that: 

1) Under the conditions precised below, a generalized 
permutator between physical states is equivalent to an 
ordinary permutator between the so called Sutherland's 
species (SS). The latter do not need to be identified with 
the physical species (PS) of the states. Indeed, the num- 
ber of PS is determined by the nature of the problem (in 
our cases they are always the four e a ), while the number 
of SS is determined by the structure of the GP entering 
the Hamiltonian. In particular, it may happen that dif- 
ferent PS constitute a single SS, so that the number of 
the latter is < 4, leading to a reduction of the dimension- 
ality of the effective Hilbert space. Through a suitable 
identification of the Sutherland's species, the first term 
of (j^) can be rewritten in the form: 

H = - ^2paN A a - oab'&ab , (4) 

A A>B 

where pa — ±1 determines the nature of the A-th species, 
even (£)/odd (O) for +1/ — 1; Naa is the number 
of neighboring sites occupied by the same species A, 
and $ab permutes objects of species A and B that 
occupy two neighboring sites, otherwise it gives zero. 
The (Jab are signs. For a given GP, the SS are to be 
identified through the subsets <Sj's of A d . In practice, 
the reduction to Sutherland's species is possible if: a) 
9 d ap = Pl Va, 13 g S t and b) 9° a/3 = o t] , Va G S, , V/3 G S, , 
where i ^ j. 

2) In the case where (|^) has <jab=+^ VA, B G £, and 
(tab independent of B for B G O and A £ £ (Suther- 
land's Hamiltonian), it is possible to reduce the number 
of even species down to only one. Indeed in this case 
the ground state energy of hamiltonian (|J) for a system 
with x even species and y odd species is equal to that 
of the same hamiltonian acting on a system with the 
same number of odd species but just one even species 
collecting all the previous ones (as implied by a simple 
extension of Sutherland's theorem, see Hj). 

Remark: For a given GP, the fulfillment of the conditions 
given at points 1) and 2) depends also on the normal- 
ization chosen to define the basis vectors. It is worth 
emphasizing that some GP, though apparently violating 
the above requirements, can be brought to fulfill them 
through a mere redefinition of the phase of a given phys- 
ical species a, i.e. \e a )j — > (— 1)^ |e a )j ■. We shall make 
use of this remark in the following. 

To illustrate how the above observations can be ex- 
ploited, we start with a known case, the AAS model [|, 
which differs from the ordinary Hubbard model only for 
a correlated-hopping term (X = 1; X = V = W = Y = 
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P = Q = 0). This model is of the form (EJ), with {? = -4, 
p, = 2, and c = — 1. Its GP reduces to (||) by identifying 
the following two Sutherland's species: A = {| f), I I)} 
(odd) and B = {|0}, | J,}}} (even). In this formalism the 
model is nothing but a free spinless fermion model, and 
its energy per site is given by 

2 

e = 2ua — 1 sin(7rn J 4) + (U — U)n-r\ — pn — c , (5) 

where ua = + n±. The phase diagram as a function 
of the filling n = N/L and U can be easily derived by 
exploiting the identity n-j-j = (n — ha)/?* and minimiz- 
ing e with respect to tia (tia £ [0, n] for < n < 1 
and n,4 € [0, 2 — n] for 1 < n < 2) and coincides with 
that derived in We also notice that the model with 
X = 1; X = 2; V = W = Y = P = Q = has the same 
energy as AAS. In fact, using the above Remark and re- 
defining the basis vector | J,|)j — » ( — l)- 7 1 |T)j, it can be 
cast in the form (||) with the same GP as AAS. 

The method just outlined allows the solution of 
a wide class of models whose Hamiltonian has the 
form (||) ||. In particular, as the AAS model dis- 
plays both a BC and a FC transitions, here we ap- 
ply it to the study of similar, but new, models in 
which further terms are included, yielding a change in 
the value of the parameters characterizing the transi- 
tion. We consider the model with coupling constants 
X = 1;X = (1 - a);Y = -a; P = -1; Q = 2 where 
a = ±1. The resulting Hamiltonian has the form (||) 
with U — —2, p = 2, c = — 1 in both cases a = ±1. 
II has diagonal entries characterized by the subsets 
Si = {1,2} , S 2 = {3} , S 3 = {4}, and off-diagonal 
entries 9°p = +1 Va G S x , (3 G S 2 and 6° af3 = a Va G Si 
or 1S2 , (3 G S3. Both the conditions a) and b) to identify 
Sutherland's species are thus fulfilled, and the species 
read: A = {\ f), | J.)} (which is 'odd' because 0*p = -1 
if a,/? G <Si) ; B = |0) ('even' because 0g 3 = +1) and 
C = I It) ('even' because 6f 4 = +1). For the case a = +1 
one can straightforwardly apply Sutherland's theorem 
(see point 2)) to reduce the number of even species to 
1, ending up with a free spinless fermion problem, where 
occupied sites are represented by A and empty sites by 
B and C. The ground state energy e per site has the 
same form as (^|), in which again njj = (n — tia)/2, and 
it has to be minimized with respect to tia- For the case 
(7 = — 1, before using Sutherland's theorem, one has to 
apply again the Remark, changing | XT) j ~* (~ I)" 7 ! IT)j- 
The expression of e is identical. 

The phase diagram is given in figj^; the lower region 
I is characterized by tla = 0, so that only doubly oc- 
cupied or empty sites are present in the ground state; 
in this region the ground state l^o) is made of the so- 
called (pure) eta-pairs, i.e. |\&o) = (?7^) |0), where 

= J2j eVp * c tf c ji = Sfe c U°l-ki witn P air momentum 
tp = 0, 7T. In the case a = +1 we have 0-pairs, whereas 
if a = —\ the pairs have tt- momentum. The latter 
case is particularly important because the 7r-pairs (and 



not 0-pairs) are expected to survive as the constraint 
X = 1 is relaxed (see ||]). In region II, delimited by 
U n _ IU = 2 — 4cos(7rn), we have the simultaneous presence 
of empty (o), singly occupied (©) and doubly occupied 
(♦) sites; this is called mixed region and the ground state 
is l^o) = (v}p) N ^ 2 \U — 00}, where \U = 00) are the eigen- 
states of the U = 00 Hubbard model Q . In both region 
I and II the ground state is superconducting, because 
the 2-particle reduced density matrix exhibits long range 
correlation §j, i.e. g(i,j) = (ty \4i c ii c ji c ji l*o> -» 
for \i — j\ — > +00. Finally, the region Ill-a (0 < n < 1) 
is made of singly occupied and empty sites; in this re- 
gion the ground-state of the U = 00 Hubbard model is 
eigenstate of the Hamiltonian and is metallic. The re- 
gion Ill-b (1 < n < 2) is the particle-hole transformed of 
Ill-a, and the metallic carriers are holes. One can show 
that at half-filling the system is an insulator with gap 
A = 17-6. 

With respect to the AAS model we observe that the 
pair-hopping term has two main effects: first it removes 
the degeneracy in (p in this region (only <p = or tt 
survive, according to the sign a of Y); secondly it raises 
the borderline of such region upwards: in fact it can 
be generally shown |J that a pair hopping term acts 
as an effective attraction (oc — renormalizing the 
Coulomb repulsion U. The superconducting region II of 
our model is enhanced also with respect to that of the 
EKS model. Indeed, although the pair-hopping term is 
also present in the EKS model (the borderlines of region 
I coincide), its effect is strongly reduced near half- filling 
due to the Coulomb attraction term between neighbor- 
ing sites (V = —1), which is known to compete with the 
formation of on-site pairs. 

As a consequence, the BC insulator-superconductor tran- 
sition occurring at half-filling corresponds to the maxi- 
mum critical value [/™ a:r = 6, higher than for all other 
exactly solved models. This is important because higher 
values of U c reduce the probability that thermal fluctu- 
ations may destroy the superconducting phase. 

Because of the particle-hole symmetry of the models 
we have considered so far, the insulating phase can exist 
just at half filling. In order to investigate FC transitions 
between finite metal-insulator regions of the phase dia- 
gram, we now discuss a simple model not particle-hole 
invariant, describing a competition between the U = 00 
Hubbard model (excluding doubly-occupancy) , and the 
pair-hopping (favouring the formation of pairs), modu- 
lated by the band parameter U (explicitly X = X = 1; 
Y = a; V — W = P = Q = 0). It is easy to realize that 
(up to the application of the Remark) this model can be 
set in the form (g) (U = -2, p = 2, and c = -1). The GP 
is now equivalent to an ordinary permutation between the 
two Sutherland's species: A = {\ j), | J.), | !|)} (which is 
'odd' because B d a0 = -1 if a, (3 G S l = {1,2,4}); B = |0) 

('even' because #33 = +1). The ground state energy per 
site is still given by eq. (||), where now [ = n — tia, ■ 
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The phase diagram -obtained by minimizing e at fixed n 
with respect to ua in the range n/2 < ua < min(n, 1)— 
is presented in fig. ^ and exhibits again four regions. 
However, due to the absence of particle-hole invariance, 
the shape is not symmetric around half-filling. 
In region I (just doubly occupied and empty sites) only 
the Y and U terms act: the model behaves like a spin 
isotropic XX model (5+ = cj^ct^, £?r = c^c^) with 
the U term acting as a magnetic field; it is well known 
that at U = Q the correlation function has a power 
law decay g(i,j) oc \i — j'l -1 / 2 , whereas g is not known 
for non- vanishing magnetic field. However, as far as 
U < 2, long-range order arises for any non-zero value 
of anisotropy. The borderline of this region is given by 
U uu = —2 cos(7rn/2). Notice that this region raises up to 
positive values of U for 1 < n < 2. The mixed region II is 
entered as the double occupancy begins to decrease from 
its maximum value, yielding the increase of the local mag- 
netic moment M = 3/4 i" 1 £V .(* |(«j,T _ "j,J.) 2 |*o) = 
3/4 (27T -1 arccos(— {7/2) — n). The value of reaches its 
minimum for ?7 n -iiia = — 2cos(7ro) when u < 1, and for 
^ii-nib = 2 when n > 1. Correspondingly, regions Ill-a 
and Ill-b are entered. The former is metallic, the ground- 
state is that of the U = oo Hubbard model, and the 
system behaves like a Tomonaga-Luttinger liquid. The 
most interesting feature is that region Ill-b is a finite 
insulating region. More precisely, at exactly half filling 
the gap is A = U — 2, while for 1 < n < 2 no empty 
site is present, and the model behaves like the Hubbard 
model in the atomic limit. Hence here the FC transition 
takes place between two finite regions, in analogy with ex- 
perimental observations on chain hole-doped compounds. 
Interestingly, for the the special value U = 2, our model 
and the L/-supersymmetric model coincide. As a conse- 
quence, our ground state energy in this case is equal to 
that obtained in pc[ . 

In this letter we have presented exact ground state 
phase diagrams of two electron models, and studied their 
BC and FC transitions. Our analysis supports the rel- 
evance of the pair-hopping term in raising the critical 
value of U for BC superconducting-insulator transitions, 
as well as the importance of particle-hole not invariant 
terms in the apperance of a finite insulating region. The 
method we used can be implemented on all those models 
described by Hamiltonian (0) in which the GP verifies 
conditions a) and 6). We stress that such GPs all corre- 
spond to integrable models [||, i.e. they are solutions of 
the Yang-Baxter Equation (consistency equation for fac- 
torizability) . The Hamiltonians exhibit therefore a set of 
conserved quantities mutually commuting. 
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FIG. 1. Ground state phase diagram of the model 
X = 1;X = (1 - a);Y = -a\P = -1;Q = 2. The model 
exhibits a BC transition insulator-superconductor transition 
at n — 1, for U c = 6. The dashed line is the EKS model, and 
[1] M.Dumm et al, Phys. Rev. B 61, 511 (2000); R. Far- the dotted line corresponds to the AAS model. 
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FIG. 2. Ground state phase diagram of the model 
X = X = 1; Y = a; V = W = P = Q = 0. AFC 
metal-insulator transition takes place for U > 2 between two 
finite regions, Ill-a and Ill-b. 
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